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The  fuel-cell  powered  bus  is  becoming  the  favored  choice  for  electric  vehicles  because  of  its  extended  driv¬ 
ing  range,  zero  emissions,  and  high  energy  conversion  efficiency  when  compared  with  battery-operated 
electric  vehicles.  In  China,  a  demonstration  program  for  the  fuel  cell  bus  fleet  operated  at  the  Beijing 
Olympics  in  2008  and  the  Shanghai  Expo  in  2010.  It  is  necessary  to  develop  comprehensive  proton 
exchange  membrane  fuel  cell  (PEMFC)  diagnostic  tools  to  increase  the  reliability  of  these  systems.  It 
is  especially  critical  for  fuel-cell  city  buses  serving  large  numbers  of  passengers  using  public  transporta¬ 
tion.  This  paper  presents  a  diagnostic  analysis  and  implementation  study  based  on  the  signed  directed 
graph  (SDG)  method  for  the  fuel-cell  system.  This  diagnostic  system  was  successfully  implemented  in 
the  fuel-cell  bus  fleet  at  the  Shanghai  Expo  in  2010. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  proton  exchange  membrane  fuel-cell  (PEMFC)  hybrid  bus 
is  completely  clean  and  environmentally  friendly.  It  is  regarded  as 
an  important  alternative  energy  vehicle  by  the  automotive  indus¬ 
try  and  several  governments.  As  a  renewable  energy  resource, 
hydrogen  can  be  produced  from  different  kinds  of  energy  sources, 
such  as  solar  energy,  bio-energy,  electrical  energy  and  other  types 
of  energy  [1].  In  July  2008  Beijing  initiated  a  clean  public  trans¬ 
portation  plan  for  the  29th  Olympic  Games.  Fuel  cell  buses  had 
operated  on  a  specific  bus  line  for  one  year  as  a  green  energy 
vehicle  demonstration  project  for  public  transportation.  For  the 
2010  Shanghai  Expo,  three  fuel  cell  buses  were  brought  into  ser¬ 
vice  to  carry  passengers  between  exposition  halls.  This  clean, 
emission-free  and  low-noise  bus  was  well  received  by  the  vis¬ 
itors  from  around  the  world.  The  State  Key  Lab  of  Automotive 
Safety  and  Energy  at  Tsinghua  University  started  research  on 
the  fuel  cell  vehicle  in  2002  and  focused  mainly  on  commercial 
vehicles  [2].  The  latest  fuel  cell  bus  fleet  provided  by  Tsinghua 


Abbreviations:  AC,  alternating  current;  AFTCU,  air  flow  and  thermal  control 
unit;  AHCU,  air  humidification  control  unit;  CAN,  controller  area  network;  CVM, 
cell  voltage  monitor;  DC,  direct  current;  DC/DC,  direct  current  converter;  FCU,  fuel 
cell  control  unit;  GSSEDM,  generalized  steady-state  electro-chemical  degradation 
model;  GSSEM,  generalized  steady-state  electro-chemical  model;  HSCU,  hydrogen 
supply  control  unit;  PEMFC,  proton  exchange  membrane  fuel  cell;  SDG,  sign  directed 
graph;  SSEM,  steady-state  electro-chemical  model;  VCU,  vehicle  control  unit. 
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University  served  at  the  2010  Shanghai  Expo  as  VIP  reception 
buses  (Fig.  1).  The  bus  is  powered  by  an  AC  motor  with  rated 
power  of  150kW  plus  75  kW  of  braking  energy  regeneration.  Two 
fuel  cell  stacks  with  a  total  rated  power  of  80  kW  and  a  80  Ah 
power  battery  are  provided  as  the  hybrid  power  source.  The 
fuel  cell  and  battery  are  separated  by  a  DC-DC  converter.  Detail 
parameters  for  the  vehicle  and  the  power  train  are  presented  in 
Table  1. 

Safety  has  been  highlighted  with  the  development  of  the  mod¬ 
ern  automotive  industry  and  has  led  to  the  mushrooming  growth 
of  fault  diagnosis  technology.  The  fault  and  incident  rates  in  mod¬ 
ern  automotive  systems,  such  as  the  fuel  cell  system,  could  result  in 
increased  casualties  and  loss  of  income,  so  fault  diagnosis  technol¬ 
ogy  has  been  playing  an  increasingly  important  role  in  ensuring 
personal  and  equipment  safety.  From  its  beginning  using  basic 
experience,  fault  diagnosis  technology  has  become  an  expert  field 
relying  on  intelligent  and  integrated  monitoring  equipment  and 
information  technology.  This  has  merged  a  large  amount  of  elec¬ 
tronic  and  computer  technology.  The  importance  of  online  fault 
diagnosis  for  the  fuel  cell  hybrid  electric  bus  is  self-explanatory 
because  of  the  high  passenger  numbers  and  long  running  times, 
as  demonstrated  over  a  year  on  a  public  transport  route.  The  fuel 
cell  system  is  a  crucial  part  of  the  overall  fuel  cell  hybrid  system, 
because  its  satisfactory  operation  influences  not  only  the  safety  of 
the  whole  bus,  but  also  its  cruising  range.  Hence,  proper  fault  diag¬ 
nosis  of  the  fuel  cell  system  during  normal  operating  processes  is 
necessary,  because  it  not  only  ensures  safety  when  the  bus  is  run¬ 
ning,  but  also  protects  the  fuel  cell  system  from  severe  damage 
under  fault  conditions. 
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Nomenclature 

Aco 

heat-exchanger  area  of  coolant  (cm2) 

Afc 

activation  area  of  fuel  cell  (cm2) 

Ca 

specific  heat  capacity  of  air  (J  kg-1  K-1 ) 

Cfc 

average  specific  heat  capacity  of  the  stack 
(Jkg-’K-1) 

Cpa 

specific  heat  at  constant  pressure  of  air  (J  kg-1  K-1 ) 

Cpw 

specific  heat  at  constant  pressure  of  the  enthalpy 
wheel  (Jkg-1  K-1) 

Qa 

average  specific  heat  capacity  of  radiator 
(Jkg-’K-1) 

Dw 

humidity  ratio  of  air  on  the  adsorbent  surface 

Eo 

open  circuit  voltage  (V) 

F 

Faraday  constant  (C  mol-1 ) 

Fv 

core  area  and  volumetric  ratio  of  the  enthalpy  wheel 

Fs 

cross-sectional  area  ratio  between  the  air  flow  and 
channel 

hCo 

convective  heat  transfer  coefficient  of  coolant 
(W  m-2 1<-1 ) 

A/lyap 

enthalpy  of  formation  when  the  reaction  product  is 
vapor  (J  mol-1 ) 

he 

fuel  cell  current  output  (A) 

i 

electric  current  in  fuel  cell  empirical  model  (A) 

Uc 

current  density  output  of  fuel  cell  (Am-2) 

Jcp 

rotational  inertia  of  compressor  (kgm2) 

I< Y 

mass  transfer  parameter 

MCm 

rotational  torque  of  air  compressor  (N  m) 

Mcp 

load  torque  of  air  compressor  (N  m) 

m 

air  flow  per  infinitesimal  unit  of  cross  section  (kg) 

mfc 

total  mass  of  the  stack  (kg) 

mNl 

nitrogen  mass  of  air  intake  (kg) 

m02 

oxygen  mass  of  air  intake  (kg) 

rcira 

mass  of  radiator  (kg) 

Nfc 

number  of  stack  cell 

Ncvm 

channels  of  cell  monitoring  system 

Pair 

manifold  pressure  (kPa) 

1  co, in 

deionized  coolant  pressure  (kPa) 

Ph 

hydrogen  pressure  (kPa) 

Qco 

coolant  heat  transfer  power  (W) 

Qfc 

thermal  power  of  fuel  cell  (W) 

Pi 

cooling  fan  relay  status  (0-1 ) 

Pu 

main  relay  status  (0-1 ) 

Ph 

hydrogen  purge  relay  status  (0-1) 

T 

stack  operating  temperature  (°C) 

Pad, in 

air  temperature  of  enthalpy  wheel  inlet  (°C) 

Pad, out 

air  temperature  of  enthalpy  wheel  outlet  (°C) 

Pah, in 

humidified  air  temperature  of  stack  inlet  (°C) 

Pah, out 

air  temperature  of  stack  outlet  (°C) 

Pa, in 

inlet  air  temperatures  of  radiator  (°C) 

Pa, out 

outlet  air  temperatures  of  radiator  (°C) 

T 

1  co, in 

stack  inlet  temperature  of  deionized  coolant  (°C) 

Pco.out 

stack  outlet  temperature  of  deionized  coolant  (°C) 

Pfc 

average  temperature  inside  the  stack  (°C) 

Pra 

surface  temperature  of  radiator  (°C) 

Pra,in 

coolant  input  temperatures  of  radiator  (°C) 

Pra, out 

coolant  output  temperatures  of  radiator  (°C) 

Vfc 

total  output  voltage  of  the  fuel  cell  system  (V) 

Ufco 

main  relay  input  voltage  (V) 

t/skl 

voltage  of  the  first  sub-module  of  the  fuel  cell  stack 

(V) 

^sl<2 

voltage  of  the  second  sub-module  of  the  fuel  cell 
stack  (V) 

Usk3 

voltage  of  the  third  sub-module  of  the  fuel  cell  stack 

(V) 

^sk4 

voltage  of  the  fourth  sub-module  of  the  fuel  cell 
stack  (V) 

^cvm,i 

cell  voltage  of  the  i  channel  (V) 

Vfc 

cell  voltage  (V) 

Wa 

air  flow  (kgs-1) 

Wco 

coolant  flow  (kgs-1) 

^N2,in 

nitrogen  flow  of  air  intake  (kg  s-1 ) 

k^N2,out 

nitrogen  flow  of  air  exhaust  (kg  s-1 ) 

Wo2,in 

oxygen  flow  of  air  intake  (kg  s-1 ) 

Wo2,out 

remaining  oxygen  flow  of  air  exhaust  (kg  s-1 ) 

Wo2,rec 

reacted  oxygen  flow  (kg  s-1 ) 

2.  Diagnosis  combining  quantitative  with  qualitative 
analysis 

There  are  two  ways  of  applying  fault  diagnosis  technology:  one 
is  based  on  modeling  and  the  other  is  based  on  experimental  data 
[3].  The  method  based  on  a  modeling  requires  a  systematic  per¬ 
formance  model  to  be  built,  and  model  changes  to  be  analyzed, 
followed  by  fault  testing  and  identification.  The  method  based  on 
data  requires  mathematical  analysis  of  a  large  amount  of  historical 
process  data  to  mine  the  information  about  the  incident.  However, 
such  a  system  is  already  in  industrial  use,  especially  in  some  process 
industries.  It  is  difficult  to  build  precise  quantitative  models,  but  the 
mechanism  influenced  by  processes  and  the  relationships  between 
elements  in  the  system  are  clear.  As  a  result,  a  new  method  based  on 
qualitative  model  fault  diagnosis  is  proposed  using  the  mechanisms 
and  relationships  to  build  a  qualitative  model  analysis  system. 

In  the  actual  situation  of  fuel  cell  control  and  diagnosis,  a 
complicated  problem  is  the  control  system  platform,  such  as  micro¬ 
processors,  cannot  afford  the  large  calculation  when  the  vehicle 
needs  a  real  time  control  to  ensure  the  vehicle  operation  and  safety, 
if  the  control  algorithm  is  an  all  model-based  process.  The  mod¬ 
els  of  the  fuel  cell  system  and  its  components  are  mathematically 
too  complicated  to  calculate  in  real  time,  especially  the  coupling 
of  different  models  or  multi-state  description  modeling.  Therefore, 
the  solution  is  using  the  large  test  data  and  empirical  formulas  to 
replace  the  large  calculation  in  real  time.  The  huge  test  data  are 
organized  as  a  lot  of  MAP  tables,  which  are  stored  into  the  non¬ 
volatile  memory  of  the  fuel  cell  control  system.  The  combination  of 
qualitative  reasoning  with  some  specific  techniques  creates  good 


Fig.  1.  A  fuel  cell  city  bus  used  at  the  2010  Shanghai  Expo. 
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Table  1 

Parameters  of  the  fuel  cell  bus  and  the  power  train. 


Parameter  (unit) 

Value 

Vehicle  mass  (kg) 

1.4672  xlO4 

Frontal  area  (m2) 

7.5 

Drag  coefficient 

0.7 

Rolling  resistance  coefficient 

oo 

X 

o 

Mechanical  efficiency  (%) 

95 

Mass  factor 

1.1 

PEM  fuel  cell  rated  power  (kW) 

80 

DC/DC  rated  power  (kW) 

80 

Ni-MH  battery  rated  capacity  (Ah) 

80 

Electric  motor  peak  power  (kW) 

150 

Electric  motor  peak  torque  (N  m) 

1.121  x  103 

Electric  motor  rated  power  (kW) 

100 

Electric  motor  maximal  rotational  speed  (rpm) 

6x  103 

opportunities  for  development,  especially  when  using  qualitative 
reasoning.  The  SDG  (signed  directed  graph)  model  uses  reasoning 
based  on  a  qualitative  model,  but  its  reasoning  methods  can  be 
divided  into  qualitative  reasoning  and  semi-quantitative  reasoning 
[4].  The  main  information  used  in  qualitative  reasoning  is  the  SDG 
model  and  its  samples,  while  semi-quantitative  reasoning  (or  semi- 
qualitative  reasoning)  relies  on  further  reasoning  based  on  the  SDG 
model  and  its  samples  by  using  quantitative  information  provided 
by  the  system.  The  processes  and  consequence  of  this  tend  to  be  a 
probability  description  or  approximate  description  to  express  the 
possibility.  The  actual  SDG  model  built  was  a  qualitative,  linear  and 
initial  responding  model,  while  a  lot  of  quantitative,  nonlinear  and 
transient  information  was  not  used  when  the  model  was  built.  If  the 
information  on  model  node  and  orientation  is  used,  it  will  approach 
the  semi-quantitative  reasoning  method.  Fig.  2  presents  the  system 
structure  of  the  fuel  cell  system  following  by  models  of  each  system, 
from  which  the  corresponding  SDG  semi-qualitative  model  will  be 
developed. 

2.1.  Stack  model 

When  optimizing  and  controlling  a  fuel  cell  system,  it  is  cus¬ 
tomary  to  describe  the  polarization  characteristics  of  the  stack  by 


using  empirical  or  semi-empirical  models.  The  most  classical  semi- 
empirical  model  is  Amphlett’s  SSEM  (steady-state  electro-chemical 
model)  [5]: 

Vfc  —  +  ^7 act  +  hohm 

<  ifect^i+fcT  +  kr  ln(c52)  +  §4rin(i)  (1) 

„  *?ohm  =  -iRmt  =  -*(£ 5  +  $6T  +  &0 

In  the  formula,  F0  is  the  open  circuit  voltage,  i  is  the  electric 
current,  T  is  the  stack  operating  temperature,  is  the  oxygen  con¬ 
centration  of  catalyst  layer,  and  §i_7  is  the  fitted  parameter.  Using 
this  basis,  Mann  and  Fowler  addressed  the  GSSEM  [6]  (generalized 
steady-state  electro-chemical  model)  and  GSSEDM  [7]  (generalized 
steady-state  electro-chemical  degradation  model)  empirical  mod¬ 
els,  respectively,  by  taking  some  factors  into  consideration,  such  as 
the  ionic  conductivity  of  the  membrane  changing  with  the  change 
of  water  content  of  the  membrane,  and  the  lifespan  of  the  cell, 
amongst  others. 

2.2.  Hydrothermal  system  model 

Theoretically  speaking,  the  thermal  power  of  a  fuel  cell  can  be 
obtained  from  the  following  formula  [8]: 

Qfc  =  %4Ac  -  Vfc)  **  Nfc'fcO  .25  -  Vfc)  (2) 

In  the  formula,  qfc  is  the  thermal  power  of  the  fuel  cell,  Nfc  is 
the  number  of  each  stack  cell,  ifc  is  the  current  density  output  of 
the  fuel  cell,  Afc  is  the  activation  area  of  the  fuel  cell,  Ahvap  is  the 
enthalpy  of  formation  when  the  reaction  product  is  vapor,  F  is  the 
Faraday  constant,  Vfc  is  the  voltage  of  each  cell,  and  Ifc  is  the  current 
output  from  the  fuel  cell. 

Vfc  can  be  estimated  using  the  measured  values  from  the  moni¬ 
toring  system  for  each  fuel  cell  voltage.  Assuming  that  Ncv m  refers 
to  the  channels  of  each  monitoring  system,  and  Vcvmi  refers  to  the 
voltage  measured  by  the  monitoring  system  in  the  i  channel,  can 
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be  calculated  from  the  following  equation: 


K  = 


i 

Ncvm 


Ncvm 

i=l 


(3) 


If  a  high  degree  of  accuracy  is  not  necessary,  Ufc  can  be  simply 
calculated  with  the  output  voltage  of  the  fuel  cell  stack  given  as: 


VfC=Wc 


(4) 


The  heat  absorbed  by  the  coolant  can  be  calculated  using  the 
following  equation  [9]: 


Qc o  —  frco^co (Tfc  —  0-5(TCO  in  +  Tco,out))  (5) 

In  the  formula,  qco  is  the  heat  transfer  power  of  the  coolant, 
hc o  is  the  convective  heat  transfer  coefficient  of  the  coolant,  Aco 
is  the  heat-exchanger  area,  Tfc  is  the  average  temperature  inside 
the  fuel  cell  stack,  and  Tcoin  and  TCOlout  are  the  input  and  output 
coolant  temperatures  of  the  stack,  respectively.  When  the  fuel  cell 
heat  exchanger  surface  is  neglected,  the  heat  transfer  equation  of 
the  fuel  cell  stack  can  be  represented  as  follows: 


mfcCfc  -jjp  =  Qfc  -  Qco  (6) 

In  the  formula,  mfc  is  the  total  mass  of  the  fuel  cell  stack,  and  Cfc 
is  the  average  specific  heat  capacity  of  the  fuel  cell  stack.  Using  the 
continuity  equation  for  fluids,  we  can  conclude  that: 

{kl^co,in  —  k^co,out  =  VVra  in  —  Wra  0ut  =  Wco 

(7) 

Wa,in  =  Wa,out  =  Wa 

The  heat  transfer  equation  of  the  radiator  can  be  described  as 
follows: 

mraOa  ^  =  WcoQo(Trajn  —  Traj0ut)  —  WaCa(Ta?0ut  —  Tain)  (8) 

In  this  equation,  mra  is  the  mass  of  the  radiator,  Ca  is  the  specific 
heat  capacity  of  air,  Cra  is  the  average  specific  heat  capacity  of  the 
radiator,  Tra  is  the  surface  temperature  of  the  radiator,  Train  and 
Tra.out  are  the  input  and  output  coolant  temperatures  of  the  radi¬ 
ator,  respectively,  Wa  is  the  air  flow,  Wco  is  the  coolant  flow,  and 
Ta  in  and  Tai0Ut  are  the  air  temperatures  in  and  out  of  the  radiator, 
respectively. 

The  flow  through  the  cooling  water  pump  is  related  to  the 
rotational  speed  of  the  pump  and  the  pressure  difference  through 
the  pump.  It  can  be  shown  that  the  flow  through  the  pump  is  in 
direct  proportion  to  the  rotational  speed,  and  the  rotational  speed 
depends  on  the  control  command.  The  relationship  between  the 
flow  at  a  given  rotational  speed  and  the  pressure  difference  can  be 
shown  by  graph  MAP,  presented  in  Fig.  3. 


2.3.  Air  system  model 

The  air  system  model  mainly  includes  the  air  compressor  and 
enthalpy  wheel  humidifier  models.  The  dynamic  characteristics 
when  the  blower  drives  the  machine  can  be  presented  as  follows: 

Jcp^=Mcm-Mcp  (9) 

In  this  equation  JCp  is  the  rotational  inertia  of  the  air  compressor 
machine,  cocp  is  the  rotational  speed  of  the  air  compressor  machine, 
Mcm  is  the  rotational  torque  of  the  air  compressor,  and  Mcp  is  the 
load  torque  of  the  air  compressor.  According  to  the  laws  of  con- 


Fig.  3.  Performance  curves  of  the  cooling  water  pump. 
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Fig.  4.  Compressor  performance  curves. 


servation  of  mass  and  conservation  of  energy,  we  can  obtain  the 
following  equations: 

d  m0o 

-gp  =  Wo2,in  -  Wo2,out  -  Wo2,rec 
<  (10) 

dniNu 

WN2’in_WN2’°Ut 

In  the  equations  above,  m02  is  the  oxygen  mass  of  air  intake, 
mN2  is  the  nitrogen  mass  of  air  intake,  W02  [n  is  the  oxygen  flow 
of  air  intake,  W02,0ut  is  the  remaining  oxygen  flow  of  air  exhaust, 
Wo2,rec  is  the  reacted  oxygen  flow,  WN2  in  is  the  nitrogen  flow  of 
air  intake,  and  WN2?0Ut  is  the  nitrogen  flow  of  air  exhaust.  When 
referring  to  the  molar  mass  of  oxygen,  the  following  equation  can 
be  developed: 

Wo2,rec  =  M02^  (11) 

The  performance  characteristics  under  steady  state  conditions 
can  be  explained  by  graph  MAP  using  experimental  detection  as 
shown  in  Fig.  4. 

Fig.  5  is  an  illustration  of  the  enthalpy  wheel  humidifier’s  the¬ 
ory,  and  shows  that,  in  general,  humidity  is  almost  100%  when  air 
is  exhausted  from  the  fuel  cell  stack.  The  enthalpy  wheel  absorbs 
heat  and  water  from  the  air  exhaust  and  stores  them  on  the  sur¬ 
face,  after  the  air  goes  into  the  enthalpy  wheel  humidifier.  When 
fresh  air  enters  the  enthalpy  wheel  humidifier,  because  the  rela- 
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DC  Motor 


Fig.  5.  Schematic  diagram  of  the  enthalpy  wheel  humidifier  [11]. 


tive  humidity  of  fresh  air  is  low,  it  will  vaporize  the  water  inside 
the  enthalpy  wheel  and  be  led  away  to  humidify  the  air  intake. 
When  the  enthalpy  wheel  humidifier  works  at  slow  speed,  the 
enthalpy  wheel  will,  step  by  step,  heat  and  humidify  the  fresh  air 
and  absorb  heat  and  water  from  the  air  exhaust  of  the  fuel  cell, 
these  being  called  the  humidification  process  and  the  regenera¬ 
tive  process  of  the  enthalpy  wheel,  respectively.  The  core  part  of 
the  enthalpy  wheel  humidifier  is  a  honeycomb  ceramic  enthalpy 
wheel.  This  is  a  corrugated  medium  composition  made  of  a  special 
high  temperature-resistant  compound,  often  with  cordierite  as  the 
main  component  [10]. 

An  infinitesimal  composition  model  of  a  humidifying  enthalpy 
wheel  can  be  built  using  the  model  building  method  for  a  desiccant 
wheel  [12,13],  and  the  mathematical  model  of  the  enthalpy  wheel 
humidifier  can  then  be  obtained  (Fig.  6). 

Therefore,  the  mass  conservation  equation  is: 


< 


dD  dD  irii  dD  _ 
~dt  +  ®  '  30  +  Ws  '  dz  " 

dW  dW  _  KYFv 

dt  30  PwFyS w 


JCyFy 

PiFs 


•(Dw-D) 


(D-Dw) 


(12) 


And  the  energy  conservation  equation  is: 

3  T  3  T  m,  3  T  aFv 

dt+C0'd^  +  Ws'dZ~  AFs(Cpa+DCpv) 


dTw  9Tw  oiFy 

~d F+®  '  “30“  =  pwFy^w  •  (Cpw  +  w  .  Cps) 


(Fw  -  T ) 
(T-Tw)  (13) 


KyFvQ 


PwFv$w  •  (Cpw  +  W  •  Cps) 


•(D-Dw) 


In  the  equations  above,  i  is  the  time  coordinate,  Z  is  the  axial 
coordinate  of  the  enthalpy  wheel,  0  is  the  circumference  corner 
coordinate,  D  is  the  vapor  portion  of  dry  air  in  moist  air,  Dw  is  the 
humidity  ratio  of  air  on  the  adsorbent  surface,  m  is  the  air  flow  per 
infinitesimal  unit  of  cross  section,  Fy  is  the  core  area  and  volumet¬ 
ric  ratio  of  the  enthalpy  wheel,  Fs  is  the  cross-sectional  area  ratio 
between  the  air  flow  and  channel,  Cpa  is  the  specific  heat  at  con¬ 
stant  pressure  of  air,  Cpw  is  the  specific  heat  at  constant  pressure 
of  the  enthalpy  wheel,  W  is  the  water  absorption  of  the  enthalpy 
wheel  ceramic,  a  is  the  thermal  conductivity  between  the  enthalpy 
wheel  ceramic  core  and  air,  I<Y  is  the  mass  transfer  parameter,  pw 
is  the  adsorbent  density,  and  <5W  is  the  depth  of  the  adsorbent. 


(a)  (b) 

(a)  relationship  between  humidification  and  speed  of  enthalpy  wheel 
(b)  relationship  between  humidification  of  enthalpy  wheel  and  speed  of  compressor 


Fig.  7.  Curves  of  humidification  function  of  the  enthalpy  wheel  during  experiments. 
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Fig.  7  shows  curve  diagrams  which  describe  the  humidification 
function  of  the  enthalpy  wheel  after  experimental  detection. 

3.  Fuel  cell  system  SDG  diagnosis  model 

Based  on  the  definition  of  graph  theory  used  in  mathematics 
[14],  and  reference  [15]  we  can  obtain  a  rigorous  mathematical 
expression  of  the  SDG  model: 

G  =  (V,E,<p,x/r)  (14) 

In  the  expression,  the  set  of  nodes  V  =  [v^,v2, ...  ,vn}  represents 
the  root  of  the  system  fault,  spur  track  £={ei,  ei, . .  .,em}  represents 
the  ‘cause  influence’  relationship  between  nodes,  map  (p\E^{+, 
-}  represents  the  positive  and  negative  effects,  map  x/r :  VW{+,  0, 
-}  represents  the  state  of  the  nodes,  <p(ek),  ekeE  represents  the 
symbol  of  the  spur  track,  x//{Vi),  Vj  e  V  represents  the  symbol  of 
the  nodes.  When  the  value  of  every  node  is  defined,  the  symbol  of 
the  spur  track  (p(ek)e{+,  -}  is  given  by  the  SDG  model,  and  the 
symbol  of  the  nodes  ^(i/j)  e  {+,  0,  -}  is  given  by  the  SDG  model 
sample. 

The  process  of  building  the  SDG  model  is  to  ascertain  the  state  of 
the  nodes  and  spur  tracks,  and  the  way  of  building  the  mathemati¬ 
cal  model  is  to  transfer  quantitative  mathematics  into  a  qualitative 
SDG  model.  Therefore,  every  assembly  is  physically  modeled  and 
tested  to  get  their  qualitative  relationships,  and  the  symbol  direc¬ 
tional  model  of  every  system  assembly  is  built  by  describing  the 


causality  between  variables  with  the  help  of  a  graphical  method, 
such  as  Fig.  8.  This  is  the  most  basic  SDG  model,  showing  that  there 
is  a  connection  between  the  node  variable  A  and  the  node  variable 
B,  while  the  symbol  between  the  lines  cp  :E^{+,  -}  refers  to  the 
relationship  that  may  have  a  positive  or  negative  effect,  which  can 
then  be  confirmed  by  models  and  practical  experience. 

Fig.  9  represents  the  qualitative  SDG  model  of  the  fuel  cell  sys¬ 
tem.  In  the  graph,  Dfc  is  the  total  output  voltage  of  the  fuel  cell 
system,  Jfc  is  the  total  output  current  of  the  fuel  cell  system,  Tahin 
is  the  humidified  hot  air  temperature  when  the  fuel  cell  stack  is 
humidified  by  the  enthalpy  wheel,  Tah,out  is  the  temperature  of  the 
air  in  the  stack  after  being  humidified  by  the  enthalpy  wheel,  Tadin  is 
the  temperature  of  the  air  entering  the  humidifying  enthalpy  wheel 
from  the  exit  of  the  blower,  Tad,0ut  is  the  temperature  of  the  sweep 
gas  from  the  humidifying  enthalpy  wheel  when  discharged  into 
the  atmosphere  after  dehumidifying,  Pair  is  the  manifold  pressure, 
PH  is  the  hydrogen  pressure,  Pco  in  is  the  deionized  coolant  pres¬ 
sure,  Tco  in  is  the  inlet  temperature  of  the  deionized  coolant  (relative 
stack),  Tco.out  is  the  exit  temperature  of  the  deionized  coolant  (rel- 


+ 


>  - ► 


Fig.  9.  Fuel  cell  system  SDG  model. 
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Notation: 

Al:  Cooling  Fan  1 

A2:  Cooling  Fan  2 

A3:  Cooling  Fan  Relay  1 

A4:  Cooling  Fan  Relay  2 

A5:  Cooling  Fan  DC/DC  Convert 

A6:  Enthalpy  Wheel  Temperature  Sensor  1 

A7:  Enthalpy  Wheel  Temperature  Sensor  2 

A8:  Enthalpy  Wheel  Temperature  Sensor  3 

A9:  Enthalpy  Wheel  Temperature  Sensor  4 

A10:  Motor  of  Enthalpy  Wheel 

All:  Hydrogen  Purge  Valve 

A12:  Hydrogen  Inlet  Valve 

A13:  Water  Injection  Valve  1 

A 14:  Water  Injection  Valve  2 

A 15:  Air  Compressor  Controller 

A 16:  Air  Compressor  Motor 

A17:  Coolant  Pump  Controller 

A 18:  Coolant  Pump  Motor 

A19:  Current  Sensor 

A20:  Protection  Diode 

A21:  Voltage  Sensor  1 


III 

L  433  | 

III 

1 

A23:  Voltage  Sensor  2 

r . 1 

t 

A24:  Power  Protection  Diode 
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A26| 
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CVM3 


CAN  B 


CVM4 


A25:  Sub  Stack  Voltage  Sensor  1 

A26:  Sub  Stack  Voltage  Sensor  2 

A27:  Sub  Stack  Voltage  Sensor  3 

A28:  Sub  Stack  Voltage  Sensor  4 

A29:  Coolant  Temperature  Sensor 

A30:  Stack  Outlet  Coolant  Temperature  Sensor 

A31:  Coolant  Pressure  Sensor 

A32:  Main  DC/DC  Converter 

A33:  Sub  Stack  Module 


Abbreviations: 

VCU:  Vehicle  Control  Unit  FCU:  Fuel  Cell  Control  Unit  CVM:  Cell  Voltage  Monitor 
AFTCU:  Air  Flow  and  Thermal  Control  Unit  AHCU:  Air  Humidification  Control  Unit 
HSCU:  Hydrogen  Supply  Control  Unit  CAN  A:  Distributed  Control  Networks 

CANS:  Cell  Voltage  Monitoring  Networks  CAN  C:  Vehicle  Networks 


Fig.  10.  Fuel  cell  control  system  and  sensor  networks. 


ative  stack),  Uskl  is  the  voltage  of  the  first  sub-module  of  the  fuel 
cell  stack,  Usk2  is  the  voltage  of  the  second  sub-module  of  the  fuel 
cell  stack,  lfsk3  is  the  voltage  of  the  third  sub-module  of  the  fuel  cell 
stack,  t/sk4  is  the  voltage  of  the  fourth  sub-module  of  the  fuel  cell 
stack,  Ufco  is  the  input  voltage  of  the  main  relay,  R\j  is  the  main  relay, 

is  the  cooling  fan  relay,  and  is  the  hydrogen  purge  relay.  The 
fuel  cell  control  system  and  sensor  networks,  which  acquire  and 
monitor  the  sensor  electronic  signals  to  provide  data  flow  for  SDG 
model,  are  shown  in  Fig.  10. 

According  to  the  definition  of  SDG,  there  are  three  states  in  every 
state  node:  xj/ :  V->{+,  0,  -},  where  “+”  means  that  the  measured 
value  of  the  state  node  is  greater  than  the  anticipated  value,  and 
similarly,  means  smaller,  and  “0”  means  it  is  in  the  normal 
range.  Valves  R i  ,R\j,Ru  represent  the  operational  state  of  the  valve; 
where  “+”  means  the  state  node  cannot  transfer  a  positive  alteration 
of  the  measure  value  of  the  node’s  correlativity,  and  means  the 
node  cannot  transfer  a  negative  alteration.  Qualitative  diagnosis 
based  on  rules  can  detect  the  occurrence  of  a  system  fault  quickly 


Fig.  11.  The  fault  transmittal  process  of  the  cooling  fan  relay. 


to  find  the  origin  of  the  fault  which  can  be  obtained  from  the  fault 
transmittal  process  in  the  SDG  model. 

Fig.  1 1  represents  a  typical  fault  transmission  process.  Because 
of  the  damage  to  the  cooling  fan  relay,  the  water  temperature  of  the 
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Table  2 

Fault  transmittal  process  of  the  fuel  cell  system  SDG  method— the  pump  loses  efficacy  due  to  leakage. 


Status 

Ufco 

Ufc 

Ifc 

Pair 

Pco,in 

Ph 

Ta 

},in  Tco,out 

Tah,in 

Tad,in 

Tah.out 

Tad, out 

Pi 

Ru 

Ph 

t/skl 

Usk2 

t/sk3 

Usk4 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

- 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

- 

0 

+ 

+ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

- 

0 

+ 

+ 

+ 

0 

+ 

+ 

0 

0 

0 

0 

0 

0 

0 

5 

+ 

+ 

- 

0 

- 

0 

+ 

+ 

+ 

0 

+ 

+ 

0 

0 

0 

+ 

+ 

+ 

+ 

Table  3 

Fault  transmittal  process  of  the  fuel  cell  system  SDG  method- 

-blower  loses  efficacy  from  being  blocked. 

Status 

Ufco 

Ufc 

Ifc 

Pair 

Pco,in 

Ph 

TC( 

},in  Tco, out 

Tah,in 

Tad,in 

Tah.out 

Tad, out 

Pi 

Ru 

Ph 

Usm 

Usk2 

Usk3 

Usk4 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

- 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

- 

- 

0 

- 

0 

0 

0 

0 

0 

+ 

0 

0 

0 

0 

0 

- 

- 

- 

- 

4 

- 

- 

- 

- 

0 

0 

0 

0 

+ 

+ 

0 

0 

0 

0 

0 

- 

- 

- 

- 

Table  4 

Fault  transmittal  process  of  the  fuel  cell  system  SDG  method- 

-hydrogen  purge  value  is  invalid. 

Status 

Ufco 

Ufc 

Ifc 

Pair 

Pco,in 

Ph 

TC( 

},in  Tco, out 

Tah,in 

Tad,in 

Tah.out 

Tad, out 

Pi 

Ru 

Ph 

Usm 

Usk2 

Usk3 

Usk4 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

+ 

0 

0 

0 

0 

2 

- 

- 

0 

0 

0 

0 

0 

0 

0 

0 

- 

0 

0 

0 

+ 

- 

- 

0 

0 

3 

- 

- 

- 

0 

0 

0 
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0 

0 

0 

- 

0 

0 

0 

+ 

- 

- 

- 

- 

fuel  cell  stack  inlet  rises  abnormally,  which  may  increase  the  water 
temperature  near  the  exit.  With  the  rising  temperature,  the  activity 
of  fuel  cell  stack  increases  and  the  voltage  rises,  so  the  output  cur¬ 
rent  flow  decreases,  as  does  the  output  power.  However,  because 
the  fault  can  be  isolated  at  the  start,  it  seldom  reaches  such  a  state. 
Usually,  when  Tco  in  exceeds  the  threshold  value,  the  control  sys¬ 
tem  will  detect  the  abnormal  state  of  the  system  and  initiate  safety 
measures  or  stop  the  engine,  and  then  confirm  or  estimate  the  fault 
cause  by  comparing  the  fault  transmittal  processes.  However,  the 
comparison  does  not  occur  in  real  time  when  the  bus  is  running, 
as  in  most  cases,  so  the  fault  will  be  alarmed  and  the  bus  will  stop 
running  to  guarantee  the  safety  of  the  fuel  cell  system  and  allow 
the  fault  diagnosis  process  to  finish. 

Tables  2-4  show  the  fault  transmittal  processes  for  different 
SDG  models. 

4.  Conclusion 

To  achieve  the  commercialization  of  fuel  cell  vehicles,  several 
obstacles  need  to  be  overcome  before  fuel  cell  vehicles  appear  on 
the  market  at  a  commercial  price.  One  of  these  is  the  development 
of  a  powerful  diagnostic  system  to  increase  the  reliability  of  fuel 
cell  power  systems  because  of  the  potential  dangers  of  hydrogen. 
In  this  paper,  a  new  fault  diagnosis  methodology,  based  on  the  SDG 
method,  has  been  presented  and  implemented  in  the  fuel  cell  sys¬ 
tems  for  shuttle  buses  at  the  2010  Shanghai  Expo.  An  advantage 
of  this  new  methodology  is  that  it  does  not  require  accurate  mod¬ 
eling  and  a  powerful  computing  capability  to  provide  a  diagnosis. 
This  can  be  carried  out  by  the  in-vehicle  control  system  and  the 
fuel  cell  control  system  during  road  use.  The  SDG  based  diagnos¬ 
tic  method  and  overall  diagnosis  protocol  were  established  during 
the  public  demonstration  of  the  fuel  cell  bus  at  the  2010  Shang¬ 


hai  Expo.  Any  system  failures  and  component  degradation  could 
be  easily  detected  and  observed  using  these  methods. 
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